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RESUME. Comme complement a une etude publiee dans ce journal 
en 2005, nous presentons des calculs explicites des champs electro- 
magnetiques selon la theorie de Maxwell moyennant autant la jauge de 
Lorenz que celle de Coulomb. On peut obtenir des expressions analy- 
tiques quand la source des champs est un dipole electrique oscillant. 
Comme avant, on trouve que les champs, qui sont calcules selon des 
methodes differentes, sont contradictoires. En outre, la raison pour le 
disaccord est decouverte: Contrairement a l'expectative les integrates 
retardees ne satisfont pas les equations non homogenes d'onde qu'elles 
sont censees resoudre. 

ABSTRACT. Complementing a study which was published in this 
journal in 2005 we present explicit calculations of fields predicted by 
Maxwell's equations both in Lorenz and in Coulomb gauge. Analytic 
expressions are obtainable, when the source of the fields is an oscillat- 
ing electric dipole. As before it is found that the fields calculated by 
different methods are at variance. In addition, the reason for the dis- 
crepancies is revealed: The retarded integrals turn out not to satisfy the 
inhomogeneous wave equations which they are supposed to solve. 



P.A.C.S.: 03.50.De; 11.15.-q; 41.20.-q; 41.60.-m 
1 Introduction 

In papers by Vladimir Onoochin [1] and by the present author [2] it was 
shown that the solutions of Maxwell's equations for point charges depend 
on the chosen gauge. Whereas in Lorenz gauge (c div A + d<p/dt = 0) 
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one obtains the well known Licnard-Wiechcrt fields, a calculation in 
Coulomb gauge (div A = 0) leads to indefinite or diverging integrals for 
the electromagnetic fields. The reason for this discrepancy was not en- 
tirely clarified, but in [2] it was suspected that an internal inconsistency 
of Maxwell's equations is responsible for the ambiguity in the solutions. 
In 1971 Donald Dunn [3] pointed out that Maxwell's differential equa- 
tions yield satisfactory results for the quasi-static instantaneous fields, 
but inclusion of the displacement current resulted in certain discrepan- 
cies which were, however, not resolved in Dunn's book. 

Our analysis in [2] led us to the suspicion that the mixture of ellip- 
tic and hyperbolic equations, which constitute Maxwell's system, actu- 
ally causes the incompatibilities between instantaneous solutions derived 
from elliptic equations (Coulomb field, e.g.) and the travelling wave fields 
as derived from hyperbolic equations. A combination of both results in 
inhomogeneous wave equations such as c 2 A(j) — d 2 cf)/ ' dt 2 — f (x, t) which 
have a strange property: The time t appearing in the d'Alcmbert oper- 
ator is the same time t which describes the temporal behaviour of the 
source. The observation point and the source may be separated by many 
light-years, but in the differential equation the events of emission and 
detection are connected as if they were happening at the same time t. Of 
course, in the "retarded" solutions one dates back the time t' at which the 
source is evaluated, by the travel time of the wave: t' =t — r/c. In view 
of this disparate conception of time in the differential equation and in 
the solution it appears, however, questionable whether the production of 
waves is correctly modelled by the inhomogeneous wave equations which 
follow from Maxwell's first order system. 

Whereas in [2] we concentrated on the production of fields by sin- 
gle moving point charges, we use in this paper an oscillating dipole of 
vanishing extension as the source for the fields. This method is due to 
Hertz who calculated his famous radiation dipole field on the same basis. 
The advantage is not only that one has analytical solutions for potentials 
and fields in Lorenz gauge, but one can also use the dipole source for 
the calculation of the fields in Coulomb gauge. The pertaining improper 
integrals do not diverge, but have well defined limiting values. 

After a recapitulation of the standard method of solution in Sect. 2, 
we introduce in Sect. 3 the Hertzian solution for a point-like dipole in 
Lorenz gauge and compare it with the solution calculated in Coulomb 
gauge. As expected from the results in [2] the fields obtained by the 
two methods are indeed different. In Section 4 we analyze a special 
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inhomogeneous wave equation with an extended source. Surprisingly 
it turns out that the standard retarded solution does not satisfy the 
given wave equation. This explains then why we have found so many 
inconsistencies in [2]. The tacit assumption that Duhamel's principle 
could be used to construct a solution of the inhomogeneous equation on 
the basis of solutions of the homogeneous equation turns out to be wrong 
in case of wave equations. The method is only suitable for the diffusion 
equation with its intrinsic infinite propagation velocity that connects 
temperature, for example, instantly with the heat source at any time. 
Finally, in Section 4 we discuss some consequences of our findings. 



2 Maxwell's equations solved for point charges in Lorenz 
gauge 

In vacuo Maxwell's first order equations are formulated as follows [4]: 

div E = 4np (1) 

f) /i 

crot£ = -- (2) 

divB = (3) 
- dE 

c rot B = 47r j + — (4) 



The potential ansatz 



1 BA 

E = =--V0, B = rot.4 (5) 

c dt w 

and the Lorenz condition 

d^~l| ,6, 
yield with (1-4) two second order wave equations of the same structure: 

A<p-^ = -4„ P (7) 

. t* 1 d 2 A 4tt 
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Lienard and Wiechert have solved them for point charges adopting 



j = pv(t) , e 



p (x , t) d 3 x 



(9) 



and obtained: 



(j> (x, t) = 



R-R-v/c 



t'=t-R/c 



(10) 



A (x, t) 



cR-R-v 



t'=t-R/c 



(11) 



where the vector R denotes the distance between the observation point 
and the position of the charge at the retarded time t' = t — R/c: 



R = x-x' (t' 



(12) 



With §5§ the fields may be derived from the retarded potentials (ITU1) and 
CLU>: 



x e R 



:R? 



1 ~^ + ^ n -^ ] - 



d v 



? -R\ 2 dt> 



B(x,t) = ^xE(S,t) , \ = l-^pn 
H He 



(13) 



Hertz has investigated the particular case of an oscillating electric 
dipole with moment ed — p(t) located at the origin. If the amplitude 
is small compared to the wavelength of the emitted radiation and small 
compared to the distance r at which the fields are observed, one may 
replace (|TT|) by the simplified expression [5]: 



A(x, t) 



p{t-r/c) 



r = \J x 2 + y 2 + z 2 



(14) 



and obtain with the Lorenz condition (j6|) for the scalar potential instead 
of (TO]): 

p (t — r/c) p-r p-f 



i (x, t) = — div- 



(15) 
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Insertion into the potential ansatz ([5]) yields for the fields of Hertz's 
dipole: 

p (9 *\ P . Z(p-r)r p 3 {P' r ) r p 

h (x, t) = -— H -g — H — — -i e2 ^ 3 

^T/^n p x r p x r , . 

B (*> t ) = ^T + ^2 ( 16 ) 

It is interesting to note that both (10, 11) and (14, 15) are actually 
particular solutions of the homogeneous wave equations ((7|) and (jHJ with 
vanishing sources p and j. The Laplace operator applied on the potential 
of a point charge is locally not defined at the charge's position so that the 
infinite density of the charge cloud does not show up in the differential 
equation. 



3 Maxwell's equations solved for point-like dipoles in Coulomb 
gauge 

In Sect. 3 of Ref. [2] it was shown that a solution of (1 - 5) in Coulomb 
gauge 

divA = (17) 

was not possible, since the improper integrals to be evaluated were ei- 
ther undetermined or would even diverge. One could suspect that this 
impasse was caused by the assumption of point charges which introduce 
singularities in the fields. If one uses point-like dipoles, however, instead 
of single charges, it turns out that the integrals in question do converge. 
This way a direct comparison of the fields calculated in both gauges (J6j) 
and (|17[) can be made. 

With p7|) and equation (TTJ) reduces to the Poisson equation 

A0 c = -4ttp (18) 

for the instantaneous Coulomb potential which has for a dipole the so- 
lution ^ 

M2,*) = £^ (19) 

Equation ([5]) becomes 

aj i c _i^ = _^J + I V ^ (20) 
c z at 2 - c c ot 
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The solution may be split into two terms Ac = A\ + A 2 where the first 
one is identical with the Lorenz vector potential (IT4l) 



Ai{x,t) = (21) 

cr 

and the second one obeys the inhomogeneous wave equation 



1 d 2 A 2 _ 1 34>c _ P(t) 3 



c z at 1 c at cr- 5 cr° 
We express the solution as a retarded integral 

a 2 (5? t) = zL fff ( ~ r '/c) 3 ~ r '/ c ) • *") ?'\ dh 



47T C _ 

V ' (23) 

in analogy to the method of solution which was applied to obtain (10, 
11) from (7, 8). Let us assume that the dipole is oriented in ^-direction. 
Because of the axial symmetry of the problem it is sufficient to evaluate 
the z-component of A 2 at y = 0: 

OO 7T 277 

a 2z = JL ( *> ( d e> [ v < 3 cos e ' - p{t - r ' /c) sin9 



Aire J J J r'R 
000 

B = ^/r 2 + r' 2 - 2 r r' (cos if' sin 6' sin 9 + cos 9' cos 9) (24) 

where we have used spherical coordinates in the integrand: 

x = r sin 9 cos f, y — r sin 9 sin (p, z = rcos9. For the cc-component of 

A 2 we obtain at y — 0: 

OO 7T 27T 

A 2x = -L J dr ' I d6 ' J d ^Ht-r'/ C ) C osf'cos9'^ 2 9' ^ 
000 

Carrying out the integration over the angles we find: 

r 00 

A 2z = J ^dr' + J ^r 4 dr', g = p(t-r'/c)(2-3 S in 2 9) 



A 2x =l Zr>sin29Ht-r>/c) dr , + f3r 2 sin29 P (t r'/c)^ 
I 10 cr d J 10 cr 
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In order to carry out the remaining radial integration we choose: 

p (t - r/c) = sinu> (t - r/c) (27) 

The integrals in (f2l))) can now be evaluated yielding in Cartesian coordi- 
nates with u) = c k: 

M (x, t) = V (J ■ V^j , r= ^ 2 +y 2 + z 2 

f(r,t) = 

(6 + 2fcV -fcV) cosk(r-ct) + (fcV (tt/2 - Si (k r)) -6)coscfci 
+kr (6-fcV) sin k (r - ct) + k 5 r 5 Ci (k r) sinckt (28) 

It is obvious that (|2"g]l. (I2T1) . and dHJ) inserted into ([S]) do not yield the 
electric held (IT5|) as given by the Lienard-Wiechert solution, since the 
sine and cosine integrals do not occur in (|16|) when (|27[) is inserted. Fur- 
thermore, insertion of (|2"5|) and (|2T|) into ([5]) does not yield the magnetic 
held as given by (fT6|) . since the rotation of (f2"Tj) results already in the 
magnetic field P^|) . but the rotation of ([2"5f does not vanish. 

These conclusions were already drawn in Sect. 3 of Ref. [2], but 
they could not be based on an explicit solution like (|2"5|). At this point 
the deeper reason for the discrepancy is still not clear, but the puzzle 
will be solved in the next Section when we take a closer look at the 
inhomogeneous wave equation. 



4 Attempt to solve an inhomogeneous wave equation 

In Sect. 6 of Ref. 2 we have used a method of solution applied to ([7]) 
which resulted in a potential at variance with (jlOp . Here we analyze the 
wave equation for the magnetic field which results from elimination of 
the electric field from ([2]) and (j4|): 

- 1 d 2 B 4tt 
A£-^r = --ro tj (29) 

In this linear equation one may split the magnetic field into two contri- 
butions 

B = B Q + B x (30) 
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where the first one obeys the Poisson equation 



4tt 

AB = rot j 

c 



(31) 



Its solution enters as a source into a wave equation for the second con- 
tribution: 



ABi - 



1 d 2 B l I d 2 B 



c 2 dt 2 c 2 dt 2 
Choosing the Hertzian dipole p(t) as the source in pip we find 



B 



p(t) x 



(32) 



(33) 



When this is substituted into (|32|) . all textbooks (e.g. [4] and [5]) suggest 
as a solution the retarded volume integral: 



Bi (x, t) 



4-7T C 3 



p(t — r' /c) x r ' d 3 r' 



(34) 



Due to the axial symmetry of the problem it is sufficient to calculate the 
y-component of at y = 0: 



OO 7T 



271 



Bu 



Air c A 



, i , , , ,p(t — r'/c) cos u>' sin 2 9' , . 
dr' / d0' / V ^ R ( 35 ) 







where R is defined in (j2~4")l . Inserting (|2T)l and performing the integration 
over the angles yields: 



Bu 



3 c 3 



r' cos lu (t — r'/c) , f r cos a; (t — r'/c) , 



dr' 



(36) 



Finally, with w = c fc one obtains for the Cartesian components of B\ 

k f (r, t) p x r 



3 r 3 p 



(37) 



/(r, i) = (l + fc 2 r 2 )cosfc(r-ci)- (l + fc 3 r 3 (tt/2— Si (fc r))) cosc/cf 
+fe r sin k(r — ct) — k 3 r 3 Ci (A: r) sin c A; t 
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It is obvious that ([3"3"f and (j3"7) inserted into (I3TJ1) do not yield the 
Lienard-Wiechert magnetic field (fl~6l) . As in Sect. 3 we encounter again 
an ambiguity of solutions. 

Furthermore, the Laplace operator applied on (|37|) yields: 



. - ,. . fc 3 , . N pxr 1 <9 2 £?n (x, t — r/c) , , 
A5i f, i) = T cos fc(r-ci)^ = 3 '-2 (38) 

Substituting this into the wave equation (|32p and integrating twice over 
time we obtain: 

B (x, t-r/c)-B x (x, t) = B Q (x, t) (39) 

This relationship is certainly not satisfied by the insertion of (|37[) . Sur- 
prisingly, it turns out that the "solution" (|3"?| does not satisfy the wave 
equation (|32|) which it is supposed to solve! 

The root of the discrepancies identified in Ref. [2] is now revealed: 
It is wrong to assume that the retarded (or advanced) functions would 
solve inhomogeneous wave equations like (I3"2l . In case of the potential 
equations (7, 8) the deficiency of (10, 11) was overlooked, since the chosen 
point sources did not explicitly enter into the differential equations as 
already remarked at the end of Sect. 2. If one chooses, however, a 
spatially extended source like (|33)) the inability of the retarded integral 
(|34[) to solve the pertinent wave equation surfaces straight away. 

It should be noted that the problem encountered is not caused by 
the splitting (13TJ1) . For the formal solution of (l3"2"j) any arbitrary function 
could be inserted into the inhomogeneity on the right-hand-side, quite 
independently of the special choice (|3"5)) imposed by the constraints of 
Maxwell's equations. 



5 Discussion and conclusion 

One may wonder how the mistake as to introduce retarded solutions into 
the framework of classical electrodynamics came about. It is certainly 
not due to Maxwell himself who did not deal with inhomogeneous wave 
equations at the price of committing a mistake (see Ref. [2] and [6]). It 
appears that Duhamel's principle, namely to construct solutions of the 
inhomogeneous equations out of solutions for the homogeneous equa- 
tions, has been wrongly applied to the wave equation. Originally the 
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method was applied to the diffusion equation dT/dt = D AT + q (x, t) 
which has a built-in property: Any change in the source has an instanta- 
neous influence on the temperature in the whole space. This may be an 
unphysical assumption, but it is built into the transport equation and 
allows Duhamel's principle to be used. There is a common time param- 
eter which describes consistently changes in the heat production and the 
instant reaction of the temperature field. 

The wave equation, however, describes travelling fields which are 
disconnected from their source the more the longer they travel. It is 
therefore not reasonable to construct solutions of the inhomogeneous 
equations with solutions that have "forgotten" their source. As pointed 
out in the Introduction and in [2] it makes no sense to connect the fields 
and their - possibly already extinguished - sources at the same time. 
Equation Q39p reveals nicely this logical inconsistency, as this relationship 
constitutes an impossible connection between the field Bq at the same 
place now and in the past. 

The question remains whether there are other methods to solve 
Maxwell's equations. Our position in this respect is the same as Dunn's 
[3] . The instantaneous quasi static fields may be calculated from elliptic 
equations such as (fT5)l or (l3~Tj) . In case of open circuits one may com- 
plement pip by a contribution containing the gradient of the Coulomb 
field. The induction law (|2|) allows determining also an instantaneous 
rotational electric field from a given instantaneous magnetic field (sec 
Appendix) . 

The travelling wave fields, however, follow from hyperbolic equations 
which must be kept separate from the elliptic equations. This cannot 
be done without committing inconsistencies which became obvious al- 
ready in Maxwell's Treatise. As long as one is not prepared, however, to 
modify Maxwell's system, this consequence is inescapable. Fortunately 
it is possible to solve the homogeneous hyperbolic wave equations with 
Cauchy-type boundary conditions that contain a physics different from 
Maxwell's equations such as Ohm's law, for example. A closed theory 
comprising both the phenomena of instantaneous quasistatic fields and 
travelling wave fields is not available yet. 



Solvability . 



11 



Appendix 

As Donald Dunn [3] pointed out, the technically important fields which 
occur in condensers, magnets, transformers etc. can be calculated for a 
given charge and current distribution from the integrals: 

E c (x,t) = IIJp(x',t) ^-^3 dV (A.l) 
=* . . -1 fffdBi(x',t) x-x' ,, , , A N 




4tt 9i 

(A.3) 

These instantaneous fields fall off with 1/r 3 and appear as sources in the 
wave equations: 

1 d 2 E w 1 <9 2 ^ 



* 19% 19 2 B, 
Ai?i^ - ~^2~Q^r~ = ~^2~Qtf~ (A.5) 

The set of equations (A.l) - (A.5) is entirely equivalent to Maxwell's 
equations (1) - (4), if we substitute 

E = E C + Ei + E w , B = Bi + B w (A.6) 

The wave equations (A. 4), (A.5), however, cannot be used close to 
the sources, since the same time appears at the observation point and 
at the source. Following Maxwell himself one must drop the right-hand- 
sides of (A. 4) and (A.5) so that one has only to solve homogeneous wave 
equations under suitable initial and boundary conditions. 

The set of equations (A.l) - (A.3) remains still useful to describe 
the instantaneous fields close to the sources p andj. In practice this 
separation of the "near field" and the "far field" world works quite well, 
but from a theoretical point of view a connection between the two kinds 
of fields would be desirable. 
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